Interaction effects in non-Hermitian models of vortex physics 
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Vortex lines in superconductors in an external magnetic field slightly tilted from randomly- 
distributed parallel columnar defects can be modeled by a system of interacting bosons in a non- 
Hermitian vector potential and a random scalar potential. We develop a theory of the strongly- 
disordered non-Hermitian boson Hubbard model using the Hartree-Bogoliubov approximation and 
apply it to calculate the complex energy spectra, the vortex tilt angle and the tilt modulus of 
(1 + l)-dimensional directed flux line systems. We construct the phase diagram associated with the 
flux-liquid to Bose-glass transition and find that, close to the phase boundary, the tilted flux liquid 
phase is characterized by a band of localized excitations, with two mobility edges in its low-energy 
spectrum. 



I. INTRODUCTION 

An external magnetic field between the lower critical 
field H c i and the upper critical field H C 2 of a type-II su- 
perconductor penetrates the material in the form of flex- 
ible flux lines (or vortices), fn high-temperature cuprate 
superconductors, which are strongly type-II supercon- 
ductors, a wide variety of thermodynamic and dynamic 
phases can be formed in a system of magnetic flux lines 
by an interplay of strong thermal fluctuations, intervor- 
tex interactions, material anisotropy, and disordering and 
pinning effects generated by various kinds of defects and 
impurities. EJ In addition to the famous Abrikosov vor- 
tex lattice phase, a number of fundamentally new liquid, 
crystalline ancLglassy phases have been proposed over the 
last decadeJaO The study of these phases and the tran- 
sitions among them is a subject of active theoretical and 
experimental research in condensed matter physics. 

A convenient way of understanding the system of in- 
teracting flux lines is provided by the formal mapping 
between the classical statistical mechanics of (d + 1)- 
dimensional directed flux lines and the nonpalativistic 
quantum mechanics of <i-dimensional bosonsB'B In this 
mapping, the flux lines traversing the sample along the 
direction of the external magnetic field, H = Hz, corre- 
spond to the boson world lines propagating through the 
imaginary time dimension. The thickness of the sam- 
ple in the z direction, L z , is analogous to the inverse 
temperature j3Ti (= h/ksT), which is the size of the 
imaginary time dimension, L T . The effective line ten- 
sion ei corresponds to the boson mass m and thermal 
fluctuations proportional to ksT play the role of quan- 
tum fluctuations proportional to h. When the flux lines 
satisfy periodic boundary conditions in the z direction 
or in the thermodynamic limit L z — > oo (corresponding 
to the zero temperature limit for the fictitious bosons), 
the boson-vortex mapping becomes precise. When free 
boundary conditions are more appropriate to the flux line 



system, the-boson mapping must be corrected by finite 
size effectsE More deta ils of the boson-vortex mapping 
will be discussed in Sec. II A and summarized in Table I. 

Using the mapping between flux lines and bosons, it 
has been suggested that, in clean high-temperature su- 
perconductors, the Abrikosov vortex lattice can melt 
and be transformed to an entangled flux liquid, which 
is an analogue of the bosonic superfluid phase.H It has 
also been argued that it is possible to produce new 
glassy phases by .adding various kinds of defects to clean 
superconductors.trt] In the present work, we focus on the 
"Bose glass" phase, which can be produced, for instance, 
by artificially creating parallel columnar defects, inside 
a superconductor using heavy ion irradiation.Qtrlj These 
defects are assumed to be distributed randomly in the 
plane perpendicular to the columns. In the spirit of the 
boson- vortex mapping, the columnar defects provide a 
static random potential for bosons, if they lie parallel 
to the external magnetic field. Interacting bosons in a 
strong static random potential can be localized and form 
the insulating Bose glass phase. The quantum super- 
fluid to Bose glass transition, occurring when the effec- 
tive disorder strength is changed across a critical value 
at zero temperature_has been studied extensively dur- 
ing the last decade£3~ll3 The knowledge obtained from 
this research can be directly applied to the flux-liquid to 
Bose-glass transition in high-temperature superconduc- 
tors with parallel columnar defects. 

If the direction of the external magnetic field does not 
coincide with that of columnar defects, it is convenient to 
separate the transverse component of the field from 
the parallel component H|| (= H\\z). When H± <C H\\, 
or equivalently when the tilt angle of the external field, 
6h = tan _1 (iJj_/if||), is sufficiently small, one can show 
that the z axis corresponds to the imaginary time di- 
mension and the transverse component H±_ plaits the 
role of an imaginary vector potential for bosons. Ej The 
corresponding fictitious quantum Hamiltonian is non- 
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Hermitian and in general possesses complex energy eigen- 
values. 

Recently, there has been great interest in non- 
Hermitian quantum mechanics of noninteracting bosons 
in a constant imaginary vector potential, which we call g, 
and a random scalar potential.EJ~t!J In the boson-vortex 
analogy, g is proportional to H^. Since the Hamilto- 
nian is non-Hermitian, the energy eigenvalues can be ei- 
ther real or complex. Hatano and Nelson have argued 
that all states with complex eigenvalues are extended, 
whereas those, with real eigenvalues are either localized 
or extendedE3 It is simpler to describe the basic phe- 
nomena in the (1 + l)-dimcnsional case. A model system 
(with periodic boundary conditions) appropriate to-flux 
lines in (1 + 1) dimensions is sketched in Fig. yjEZI In 
the absence of the imaginary vector potential, it is well- 
known that all eigenstates are localized (in both one and 
two dimensions) and the corresponding eigenvalues are 
real. This continues to be the case if g is sufficiently 
weak or the random potential is sufficiently strong. As 
the imaginary vector potential, or the transverse mag- 
netic field, is increased above a critical value gi, there 
appear a few extended states with complex eigenvalues 
in the center of the energy band. The spectrum of this 
group of extended states is shaped like a bubble on a 
complex energy plane and the two points where the bub- 
ble is attached to two real lines of localized states are the 
mobility edges separating the extended states from the 
localized ones (see Fig. |^). As g increases further, more 
states change from being localized to extended and the 
complex bubble expands. Ultimately, above another crit- 
ical value <?2 of the imaginary vector potential, all states 
including the ground state become extended and the en- 
tire spectrum is shaped like a single complex bubble. The 
basic physics is similar in two dimensions, although lo- 
calized and extended states coexist near the band center 
and the spectrum is much more complicatedO 

These results suggest the existence of a novel dereal- 
ization transition even in one- and two-dimensional non- 
Hermitian systems. In a small external transverse mag- 
netic field, the flux lines remain in localized states and are 
strongly pinned by columnar defects. Since the system 
displays perfect diamagnetism against the transverse field 
in this case, this phenomenon has been termed the trans- 
verse Meissner effect. However, above a certain critical 
transverse field, the flux lines are depinned from colum- 
nar defects and form a tilted flux liquid. A schematic 
phase diagram of high-temperature superconductors in 
the T-H± plane is shown in Fig. f| 

In spite of the large number of papers devoted to the 
non-Hermitian problem, there has been very little re- 
search on the effects of interactions on the derealization 
phenomena in non-Hermitian systems. One may attempt 
to take account of the interaction effects by assuming 
that the nature of the spectrum in the interacting case 
remains the same as in the noninteracting one and forbid- 
ding multiple occupancy of the single-particle localized 
eigenstates. Let us consider the (1 + l)-dimensional case 



and assume that the transverse field H± is in a range such 
that there is a mobility edge separating low-energy local- 
ized states from high-energy extended ones. In the boson- 
vortex mapping, the parallel component of the field, H\\ , 
corresponds to the chemical potential [i of the bosonic 
system. As H» is increased at a fixed H±, the states are 
filled in order of increasing energy to obtain the ground 
state. At a critical value of Hn , all states below the mo- 
bility edge are filled and additional flux lines begin to fill 
the extended states just above the boundary. These ex- 
tended states describe tilted flux lines and, the tilt angle 
is known to be finite at the mobility edge.E-0 The Hartree- 
Bogoliubov results presented here suggest that this naive 
picture, which is similar to the "inert layer" picture used 
to explain the superfluid onset transition in-iHe, adsorbed 
in porous media or on random substrates,Ein23 is wrong 
and needs to be substantially modified. 

Hwa et al. have developed a physical picture of the de- 
localized phase and the derealization transition in (1 + 1) 
dimensions in the presence of interactions E3 When H± 
is larger than a critical value, each flux line jumps from 
one columnar defect to another via kinks. These kinks 
line up and form chains to reduce the interaction en- 
ergy cost. The density of the chains goes continuously to 
zero as H± decreases to the critical value and the univer- 
sality class of the derealization transition is predicted 
to be the same as that of the standard commensurate- 
incommensurate transition in (1 + 1) dimensions. The 
physical picture in the (2 + l)-dimensional case is less 
clear and expected to be more complicated than that 
in the (1 + l)-dimensional one. More recently, Lehrer 
and Nelson studied the non-Hermitian Mott-insulator to 
flux-liquid transition occurring in a system of interacting 
flux lines with a periodic array of columnar defects, us- 
ing both the mean-field theory and the renormalization- 
group method.E3 The generic universality class of the 
non-Hermitian transition was found to be the same as 
that of flux lines entering the Meissner phase at the lower 
critical field H c \- 

The main aim of this work is to develop a quantita- 
tive theory for understanding the delocalized phase and 
the derealization transition in interacting non-Hermitian 
systems in the presence of strong disorder. We apply 
the theory to the disordered flux liquid phase and to the 
flux-liquid to Bose-glass transition in high-temperature 
superconductors with columnar defects. Generalizing 
the model of Ref. 17 to the short-range-interacting case, 
we introduce a d-dimensional non-Hermitian Hubbard 
model for lattice bosons in a random potential: 

+ ^x4 fl x + 44 a xfl x , (i) 

X X 

where 4 and a x are boson creation and annihilation op- 
erators at lattice site x respectively. The hopping integral 
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t sets the scale of the boson "kinetic energy" and U is the 
strength of the repulsive on-site interaction. The random 
potential V x is assumed to be uncorrelated in space and 
uniformly distributed between — A and A. The constant 
vector g is a non-Hcrmitian external field and the vectors 
{e„} are the unit lattice vectors. We assume a hypercubic 
lattice with the lattice constant a, therefore, | | = a. It 
is convenient to introduce a dimensionless non-Hermitian 
field h = ga/h. The connection between this model and 
the physics of vortices will be discussed in more detail 
in Sec. II A. The transverse component of the external 
magnetic field Hi is related to the non-Hermitian field 
h via 



4>qtl±a 

4irk B T ' 



(2) 



where 0o = wTic/e is the magnetic flux quantum. We use 
boson notation in the discussion below. 

The Hermitian version of the boson Hubbard model 
(|l|), which is obtained by setting g = 0, and its contin- 
uum counterpart have been studied extensively using a 
variety of methods, which include numerical simulations, 
scaling analysis, renormalization group calewlfttipps and 
the Hartree-Bogoliubov approximation.! 1 °n 1 3'ETEj In- 
spired by an earlier work by Lee and Gunn,E3 Singh and 
Rokhsar developed a numerical method of studying the 
disordered superfluid phase and the quantum superfluid 
to Bose-glass transition-in a system of interacting bosons 
in a random potential,EEI for the purpose of understand- 
ing some experiments done with 4 He adsorbed in porous 
media.EirEj Using a combination of the Hartree approx- 
imation and the Bogoliubov theory, they calculated the 
condensate depletion and the superfluid density, as well 
as the excitation spectrum and the density of states. 

In the present work, we generalize Singh and Rokhsar's 
method to the non-Hermitian case. The first step is 
to determine the ground state wave function in a self- 
consistent manner using the Hartree approximation. We 
find that the random potential is screened by a (Hartree) 
repulsive interaction with other particles in the ground 
state. The screened potential is weaker and smoother 
than the original random potential. In the next-step, we 
apply the Bogoliubov canonical transformation^ to the 
ground and excited states obtained using the Hartree ap- 
proximation and derive the Bogoliubov excitation spec- 
trum and the ground state energy. By taking derivatives 
of the ground state energy E g with respect to the non- 
Hermitian field, we calculate the quantities such as the 
imaginary current J/ and the superfluid density n s . 

The imaginary current in the ground state is defined 

by 



d_ 
<9g 



Eg 

N 



(3) 



where N is the number of bosons (or flux lines) This 
quantity corresponds to the transverse component of the 



magnetic induction, B^, in flux line systems. The vortex 
tilt angle, 



= tan 



i (B± 
Bn 



is related to Jj by 



Jj ta 
J ti 



(4) 



(5) 



The quantity Jo can be expressed as ha /m, since, in 
the continuum limit, the hopping integral t is related to 
the boson mass m (this is the tilt modulus of the vortices) 
by 



ta 2 ' 



(6) 



In a system of volume V, the superfluid density is de- 
fined by 



h2 Hm 32 r^(g+<%) 



ta 2 s g ^o d{5gf 



V 



d 2 



sh^o d{5h) 2 



Eg (h + Sh) 

tv 



(7) 



where 5h = (a/h)5g. The direction of the vector 5g is 
that of the superflow. In general, the superfluid density 
depends on the relative angle between the non-Hcrmitian 
field g and the superflow. The definition (Lai agrees pre- 
cisely with the more conventional definitiorH of n s based 
on the dependence of E g on periodic phase variations (of 
wavevector kg) of the wave function: 



m d 
<y hm — 
% 2 k ^Q dk^ 



E g (k ) 



V 



(8) 



It turns out that the superfluid density is directly pro- 
portional to the vortex part of the inverse tilt modulus, 
c| 4 , in flux line systems: 



(9) 
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where n is the density of flux lines .E3 In three-dimensional 
superconductors, n s and n represent areal densities. 

Both J 1 and n s are zero in the Bose glass phase and 
nonzero in the superfluid (or flux liquid) phase. By com- 
puting these quantities for various parameter values, we 
are able to construct a phase diagram associated with 
the Bose-glass to flux-liquid transition. When combined 
with the results on the Hartree and Bogoliubov spectra 
and associated wave functions, these calculations provide 
a physical picture for the delocalized phase and the der- 
ealization transition substantially different from the non- 
interacing case, as discussed below. 

Let us again consider the one-dimensional problem 
where there is a mobility edge separating low-energy lo- 
calized states from high-energy extended ones. As H\\ 
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(corresponding to the boson chemical potential) is in- 
creased from the Bose glass side, the localized states are 
filled in order of increasing energy. At the same time, 
the interaction effect, which is proportional to the den- 
sity and screens the random potential, becomes stronger. 
Since the screening effect is most effective for the many- 
body ground state, this ground state becomes extended 
one at a critical value of H», well before the chemical 
potential hits the mobility edge in the single particle pic- 
ture, and the system becomes a tilted flux liquid. As H\\ 
is increased further, more low energy states are trans- 
formed to extended ones and eventually all low-energy 
states will become extended. These arguments suggest 
that, close to the Bose glass phase boundary, the tilted 
flux liquid is characterized by a low-energy band of Bo- 
goliubov excitations delocalized by hopping, a band of 
localized excitations, and then a band of more conven- 
tional delocalized states. Its spectrum thus contains two 
mobility edges. 

The outline of the paper is as follows. In Sec. [il], we 
develop the basic formalism using the Hartree and Bo- 



goliubov approximations. In Sec. Ill, we present the nu- 



merical results in one dimension on the energy spectra 
and associated eigenfunctions both in a random poten- 
tial and in the presence of a single impurity. The nature 
of the eigenfunctions is studied by computing the par- 
ticipation ratios and the winding numbers. In Sec. 
we present the numerical results on the vortex tilt an- 
gle and the tilt modulus of (1 + l)-dimensional flux line 
systems both in a random potential and in the presence 
of a single impurity. We also construct the phase di- 
agram associated with the flux-liquid to Bose-glass (or 
Bose-insulator) transition. In Sec. [V], we conclude the 
paper with some remarks. In Appendix |X[ we develop 
an analytical theory of both continuum and lattice non- 
Hermitian models of a Bose gas in a weak random po- 
tential using the Bogoliubov approximation. 



II. NON-HERMITIAN HARTREE-BOGOLIUBOV 
THEORY 

A. Model 

Our main interest is in understanding the physics 
of flux lines in high-temperature superconductors with 
columnar defects. Therefore it is useful to clarify the con- 
nection between the boson Hubbard model (Q) and the 
model free energy F for N flux lines in a sample of thick- 
ness L z in the z direction (perpendicular to the Cu02 
planes) in the presence of M columnar defects aligned 
in the z direction, but randomly distributed in the xy 
plane: 



„ N 

f=-t 



dr^z) 



dz 



dz 



]T / 'v D [n(z)]dz 



iH 



E 



* dYj(z) 

4ir ^ Jn dz 

M 

V D [r l (z)]=J2Vi[\r l (z 

k=\ 



dz, 



(10) 



where and Rfc are two-dimensional vectors in the xy 
plane. This free energy is a good approximation if the 
mass anisotropy m z /m± in the underlying Ginzburg- 
Landau model is sufficiently large, as is the case in 
cuprate superconductors.u U(r) is the repulsive interac- 
tion potential between flux lines, which can be taken to 
be local in z, while V\ is the (attractive) interaction be- 
tween a flux line and a columnar defect. The effective 
line tension 



e\ = — - ln/t, 

m z \ AttX J 



(11) 



where k = A/£ is the ratio of the London penetration 
depth A and the coherence length £, arises from the tilt 
energy of the lines. 

The canonical partition function for a system of N lines 
is given by the path integral 



Z = f[Jv[r z (z)}e^ F ^^, 



(12) 



which can be rewritten in terms of the imaginary-time 
evolution operator e - L ^ n / h as 



Z 



(13) 



where the bra and ket vectors are the initial and final 
states, respectively. The lattice version of the quantum 
Hamiltonian appearing in the above equation is precisely 
the boson Hubbard Hamiltonian (|I|). The parameters of 
this model, expressed in the original flux line language, 
are summarized in Table || 

In the main part of this paper, we will concentrate on 
the (1 + l)-dimensional case, though the theory devel- 
oped here can also be applied to (2 + 1) dimensions. For 
convenience, we rewrite the Hamiltonian (11) as 



H 



u 



+ J2(V t -^)a] ai + — Y^a\a\aia u (14) 

i i 

where h — ga/Ti. We assume that the total number of 
lattice sites is L. The density is then equal tons N/L. 
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B. Hartree approximation 

In general, right eigenfunctions of a non-Hermitian op- 
erator are not simple complex conjugates of the corre- 
sponding left eigenfunctions. Let us consider the right 
Hartree state 

i> R (rur 2 ,---,r N ) = ct>$(n)<P$(r 2 )---c(>$(r N ), (15) 

where all bosons are condensed into the same state de- 
scribed by the normalized single-particle right eigenfunc- 
tion (f>Q"(ri), as a variational ground state for (^4[). The 
corresponding left Hartree state is 

ip L {r U r 2 ,---,r N ) = $ (n) $ (r 2 ) ■ ■ ■ $ (r N ) . (16) 

The many-particle states ([Tj]) and (|l6| ) are explicitly 
symmetric under particle exchange. 

The expectation value of the Hamiltonian ( |l4| ) between 
the left and right Hartree states is 

(^|(w-/xiv)hO = 

- ^ E t e ^o ( l + + e- h $(i)<P$(i + 1)] 

i 

i 

i 

In order to minimize Eq. ( p"7| ) with respect to the single- 
particle state </>o(*)> w e need to solve the discrete nonlin- 
ear Schrodinger equation: 

-~ [e h ^(i - 1) + e- h <^(i + 1)] + WrfZii) 

= (A« + eA)0?(i), (18) 

where the effective single-particle potential Wj is defined 
by 

W i = V i + UN^{i)<f> R (i). (19) 
Similarly, we find that the wave function 4>x(i) satisfies 

~\ [e- h ti(i - 1) + e h ^(i + 1)] + Witf(i) 

= (A* + e A )^(t). (20) 

00 (*) an< i 0o are the ground state wave functions of 
Eqs. (jl^) and (jSTJ). We assume that the wave functions 
satisfy the periodic boundary conditions appropriate to 
Fig. | 

<ft(L) = <ft(0), ^(L + l) = ^(l), 

^(i) = ^(0), ^(£ + l)=^(l). (21) 

Due to the simplicity of the Hamiltonian and the 
boundary conditions, the eigenfunctions (j\ ,L {i) andr-the 
eigenvalues e\ satisfy simple symmetry relationshipsE-3 



cf> L x (i;h) = [^(i;-h)]*, 

ex(h) = [e x (-h))*. (22) 

The chemical potential fi is chosen so that the (real- 
valued) eigenvalue eo is equal to zero. This condition 
ensures that the Bogoliubov quasiparticle spectxum has 
no gap, in conformity with general theorems. lB From 
Eq. (jfj) or @, we then find 

» = -\ E [ e ^o (* + 1)0?(O + e- h <ft{i)${i + 1)] 
+ E + ^ E W 2 • (23) 

The ground state 0q : («) (0o (*)) an d the excited states 
{0A^o(*)i ({0a#o(*)D form an orthonormal basis for 
single-particle states. In Sec. |lj], we Wlu solve Eqs. ( [Ls| - 
pf]| ) numerically in a self-consistent manner. 

C. Bogoliubov approximation 

In the Bogoliubov approximation, one assumes that 
the ground state (or the condensate) is macro scopically 
occupied with occupation number Nq. This assumption 
is not true in (or sufficiently close to) the Bose glass 
phase in any dimension. Even though the superfluid to 
Bose-glass transition is possible in (1 + 1) dimensions, 
phase fluctuations will actually lead to algebraic decay 
instead of long-range order in the boson order param- 
eter. We illustrate this behavior in Fig. ||, where we 
plot the occupation number of the Hartree eigenstates, 
N\ = (b x b x ) (see Eq. (p6|)), measured in the Bogoliubov 
ground state against the winding number N w of the cor- 
responding Hartree eigenfunctions. We find that the oc- 
cupation number rises sharply as the ground state with 
N w = is approached. Though the ground state is not 
macroscopically occupied, it still plays a dominant role in 
long-wavelength physics and the Bogoliubov approxima- 
tion can be used to calculate properly quantities such as 
the imaginary current and the superfluid density.LS Our 
theory is basically a mean-field approach and is expected 
to be valid away from the critical region. 

From now on, sums over A (= 1, 2, • • • , L — 1) will al- 
ways exclude the condensate. We expand the boson field 
operator at in the complete set of right annihilation op- 
erators bff and {b^}: 

a i = $(i)bK + '£<l>x(i)bl (24) 

A 

Similarly, we expand a\ in terms of the left annihilation 
operators 6q and 

4 =0oW&o +E^«&A- (25) 

A 
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In the Bogoliubov approximation, we replace the oper- 
ators 6q and by y/No, where No is the number of 
bosons in the condensate. The total number of bosons is 
the sum of those in the condensate and those not in the 
condensate: 



N = N + Y,b L xbl 



(26) 



Expanding to first order in the depletion of the conden- 
sate, we get 



'No 



N 



2yfN ^ 



(27) 



Substituting Eqs. (|2J), (|25|) and (|27|) into Eq. and 
keeping all terms second order in fr^and b x , we obtain a 
quadratic Bogoliubov-type Hamiltonian 

i 

+ ]T (axx^, + \B xx ,b L x b L x , + icAv&X/) . (28) 

where the matrices A, B and C are defined by 

A X y = e x S xy +UNj2<l>o (WoiWx (i)<t>x> (*)> 

i 

i 

C XX , = UNj2[4(i)} 2 <t>x(i)<P*(i)- (29) 

The matrices B and C describe the processes (includ- 
ing effects due to the random potential!) where a pair 
of particles are scattered out of and into the condensate 
respectively, while the second term of the matrix A de- 
scribes single-particle scattering by the condensate. 

We diagonalize the Bogoliubov Hamiltonian ( p8| ) by 
introducing a set of quasiparticle creation and annihila- 
tion operators 7^ and 7 M (fi — 1, 2, • • • , L — 1), defined by 
the non-unitary canonical transformation 



A 



(30) 



where the coefficients u m a and i> m a arc defined by 

u^xQi) = [u^\(-h)]* , v^x{h) = [v^x(-h)]* . (31) 
The quasiparticle operators satisfy 



(32) 



where w M is the (complex- valued) quasiparticle excitation 
energy. The transformation coefficients and turn 
out to obey the eigenvalue equations 



^ {Axx'U^x- + Bxx'V^x') = u^u^x, 

A' 

^ (CaA'W^A' + AvAVpA') = -WpW^A 
A' 

and the orthonormality conditions 

^ (u^xu^X - VVa) = <W- 



(33) 



(34) 



In Sec. Ill, we solve these equations numerically to ob- 



tain the coefficients u^x and v^x and the eigenvalues 



Finally, using Eq. (23), we find that the diagonalized Bo- 
goliubov Hamiltonian has the form 



H B = E g + 5^^7+7^, 



E, 



g 
N 



[eV L (* + 1)^(<) + e- h rt(i)cj>K(i + 1)] 



C7n 



+ 5>^)0 o fl ( l ) + — L]T [^(*)^(i)] 



2n£ 



where -E^ is the ground state energy. 



(35) 



III. SPECTRA AND WAVE FUNCTIONS 

A. Ground state wave function and screening 

The first step of our numerical calculation is to solve 
Eqs. (|l8|-p0|) with the periodic boundary condition (pi]). 
In all of our calculations, we will fix the hopping integral 
t to a constant value. Then there are three natural di- 
mensionless parameters in Eqs. ( |l8| - |20| ), that is, the non- 
Hermitian field h, the interaction parameter Un/t and 
the disorder parameter A/t. It turns out that the effect 
of disorder is negligible if h or Un/t is sufficiently large, 
or of course if A/t is sufficiently small. In these cases, 
we begin with (normalized) trial condensate wave func- 
tions 4>q(i) — <f>o"(i) = 1/vL, which is uniform in space. 
Then we solve the nonlinear Schrodinger equations with 
the screened random potential Wi = Vi + Un numerically 
to find new left and right single-particle eigenstates. We 
obtain updated trial condensate wavefunctions by mix- 
ing the initial guess with the ground state eigenstate of 
the screened potential, using the so-called Broyden mix- 
ing method widely used, in self-consistent computations 
of nonlinear equations .Ell This procedure is repeated until 
a convergence is achieved. 



G 



In the presence of disorder, the ground state wave func- 
tion is nonuniform in space. We use the solution obtained 
for a large h or Un/t, or for a small A/t as trial wave 
functions for the cases with a smaller h or Un/t, or with 
a larger A/t. The convergence is very slow if the system 
size is large and the effective disorder is strong. For all 
numerical results presented here, the lattice size is lim- 
ited to L < 200. 

In Fig. ||, we plot the local condensate density n (i) — 
4>q (i)4>Q(i) in a single realization of the random potential 
for h = 0.1, Un/t = 0.5, A/t = 1 and L — 60. We notice 
that the condensate is indeed highly nonuniform. In fact, 
the local condensate density is large (small) at the local 
minima (maxima) of the bare random potential Vj as 
shown in the inset of Fig. ||, so that the screened random 
potential is smoother than Vi- In Fig. ^, we compare 
the screened random potential for the same parameter 
values as in Fig. [5] with the corresponding bare random 
potential. The screening is more effective for larger values 
of Un/t. We have verified that the screening does not 
produce long-range correlations in the random potential, 
as illustrated in Fig. [?]. 

B. Hartree and Bogoliubov spectra 

After calculating the self-consistent eigenvalues and 
eigenf unctions using the Hartree approximation, we use 
them to solve the Bogoliubov eigenvalue equations ( |33| ) 
and (|i|). The eigenvalues w M obtained from this calcu- 
lation give the Bogoliubov quasiparticle spectrum. 

In Fig. ||(a), we show the Hartree spectra of a 200- 
site lattice in a single disorder configuration for Un/t = 
0.5, A/t — 1 and h — 0.8,0.6,0.3. For these parameter 
values, the system is deep inside the superfluid region. 
When the non-Hermitian field is strong, the spectrum is 
close to that of a clean system, with no real eigenvalues 
except at the very bottom and top of the band. All states 
are extended in this case. As the non-Hermitian field h 
is reduced, some of the eigenvalues near the top of the 
band become real at some critical value of h. This band of 
real-valued eigenvalues corresponding to localized states 
grows continuously as h is further reduced. We note that 
there is a mobility edge separating low-energy extended 
states from high-energy localized states. 

In Fig. ||(b), we show the Bogoliubov excitation spectra 
for the same parameter values as in Fig. |8|(a). Similarly 
to the Hartree spectra, there is a low-energy complex 
bubble corresponding to a band of extended states. For 
small values of h, a band of (localized) states with real 
eigenvalues appears at the top of the bubble. The Bogoli- 
ubov spectrum shows a linear behavior near the ground 
state, similar to results for non-Hermitian "superfluids" 
without disordered This linear vanishing of Im ui^ with 
Re ujfj, is due to the collective phonon excitations charac- 
teristic of the superfluid phase. We notice that the effect 
of disorder is weaker in the Bogoliubov spectra than in 



the Hartree spectra in the sense that the critical value of 
h where localized states begin to appear at the top of the 
band is smaller. 

Next, we consider the change of the Hartree and Bo- 
goliubov spectra as we cross the Bose glass-supcrfluid 
phase boundary by increasing the value of h, starting 
from the Bose glass side. When h is zero, all eigen- 
values are real and, in a sufficiently disordered system, 
low-lying eigenstates including the ground state are lo- 
calized. As h is increased, a bubble of extended states 
with complex eigenvalues appears in the center of the 
band, but the low-energy states with real eigenvalues re- 
main localized. As h crosses the phase boundary, the 
ground state becomes extended and the system is trans- 
formed to a superfluid. Within the mean-field theory 
framework, the condensate is formed, even though it ex- 
hibits a slow algebraic decay at T = in (1 + 1) dimen- 
sions. As h is increased further, a few lowest-lying states 
become extended and form a small bubble in the complex 
plane. Note that there is a small band of localized excited 
states sandwiched between two bands of extended states 
in this parameter regime! The entire spectrum thus con- 
tains three mobility edges. In the continuum limit where 
lattice effects can be neglected, we expect the two low- 
energy mobility edges survive, while the top one disap- 
pears. For larger values of h, two bubbles grow in size 
and eventually merge. This "double bubble" scenario is 
illustrated in Fig. [|for a single disorder configuration and 
for Un/t = 0.05, A/t = 1, h = 0.3, 0.44, 0.5 and L = 60. 
Again, the effect of disorder is weaker in the Bogoliubov 
spectra than in the Hartree spectra, and the character- 
istic linear behavior in the Bogoliubov spectrum is ob- 
served in the small bubble including the ground state. 
We have checked that this "double bubble" scenario is 
quite robust and occurs in any realization of the random 
potential. 

C. Participation ratio and winding number 

The nature of the Hartree eigenstates can be studied 
by calculating the participation raticEJ 



of state A. This quantity is a measure of the number of 
sites where the wave function is not negligible. In case of 
plane waves, P\ is equal to the total number of sites and 
diverges in the thermodynamic limit. For more general 
extended states, a weaker condition 

lim In P A cx In L (37) 

L — >OG 

is applied. For a localized state, P\ is proportional to 
the localization volume and approaches a constant in the 
thermodynamic limit. 
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In Fig. 10, we plot the participation ratios for Un/t — 
0.5, A/t = 1, h = 0.8,0.3 and L = 100,200. We find 
that, for h — 0.8, the participation ratios for all states 
scale with the size of the system. This can be considered 
as evidence that all eigenstates in this case are extended. 
When h = 0.3, the participation ratios for the states 
having real eigenvalues at the top of the band do not 
scale with the size of the system, while those for the states 
with complex eigenvalues and the ground state do. We 
conclude that the states near the top of the band with 
real eigenvalues are localized and there exists a mobility 
edge separating them from low-energy extended states. 

A better criterion for determining the nature of 
Hartree eigenstates in one dimension is provided-by the 
winding number of the Hartree eigenfunctionHa This 
winding number is the generalization for non-Hermitian 
disordered systems of states labelled by wavevectors in 
translationally invariant problems. When the eigenvalue 
e\ is complex, the eigenfunction (f>^' L (i) rotates an inte- 
ger times around the origin on the complex plane, as the 
spatial dimension (which is a circle) is traversed once. 
This integer rotation or winding number N w is well- 
defined even in strongly-disordered cases. It has been 
argued that the states with nonzero winding numbers 
are always extended Ba 

In Fig. [n](a), we compare the Hartree spectrum and 
the winding numbers for h = 0.3, Un/t = 0.5, A/t = 1 
and L = 200. In Fig. |ll](b), we show the curves for 
h = 0.44, Un/t = 0.05, A/t = 1 and L = 60. A bubble of 
extended states including the ground state is clearly seen. 
If the eigenvalue is real, the corresponding eigenfunction 
is also real, and the winding number is undefined. In this 
case, we arbitrarily set N w = 0. 



D. Non-Hermitian vortex physics with a single 
impurity 

Some features of the Hartree spectrum can be better 
understood by examining the cases with a single impu- 
rity. In Fig. |l^(a), we show the spectra in the presence of 
a single attractive impurity with the strength Vd = — 5i. 
In the noninteracting case_there is a single bound state, 
which is the ground state.E3 As the interaction parameter 
Un/t increases, this bound state moves toward the band 
of extended states with higher real parts of the eigenval- 
ues and, at the same time, the highest energy eigenstate 
becomes a bound state. Above some critical value of 
Un/t, the ground state merges into the band of extended 
states and another bound state appears at the top of the 
band. 

This behavior can be understood from_the shape of 
the screened potential Wi shown in Fig. 
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(b). As Un/t 

increases, the potential well of the screened potential be- 
comes shallower and the ground state changes from a 
bound state to an extended one. On the other hand, 
the inverted screened potential, — Wi, has two potential 



wells, which can be sufficiently deep to support one or two 
bound states for intermediate values of Un/t. This ob- 
servation explains the appearance of hole bound states at 
the top of the band. Extending the present argument to 
the cases with many (but dilute) impurities with random 
strengths, we can qualitatively understand the Hartree 
spectra in the presence of a random potential shown in 
Figs. H and ^. 

In Fig. [13| , we show the spectra and the screened po- 
tential in the presence of a single repulsive impurity with 
the strength Vd — 5t. As Un/t increases, the hole bound 
state with the highest energy moves toward the band of 
extended states and merges with it at a critical value of 
Un/t. This can also be understood from the shape of the 
screened potential. 



IV. TILT ANGLE, TILT MODULUS AND THE 
BOSE-GLASS TO FLUX-LIQUID TRANSITION 

After obtaining the Hartree eigenvalues and eigenfunc- 
tions and the Bogoliubov eigenvalues, we use them to 
compute the ground state energy per particle, E g /N, 
given in Eq. (|a) . We note that the last term of Eq. fl35| ) 
has an explicit dependence on the density n, while the 
first three terms do not. This suggests that the density is 
another independent parameter in addition to Un/t, A/t 
and h. For all parameter values used in this work, the 
last term of Eq. (|35| ) is numerically much smaller than 
the other terms, unless the density is very small. In the 
rest of our calculations, we will fix n = 0.6. 

At the final stage of our calculations, we take numerical 
derivatives of the ground state energy with respect to the 
non-Hermitian field h and calculate the vortex tilt angle, 
9 V , and the vortex component of the inverse tilt modulus, 
c v 44 -\ using Eqs. (g), (0) and ®. In Fig. [u|, we plot 



9„ and 



44 



as functions of h, which is proportional to 



the transverse external field H±, for different values of 
Un/t. It is clearly seen that there is a transition from 
the Bose glass phase with 9 V = and c| 4 = oo to the flux 
liquid phase with both quantities finite, as h is increased 
through a critical value h c . The transition is broader 
when the interaction is stronger. It has been predicted 
in Ref. 24 that, close to the phase boundary, c v iA _1 ~ 
[h — h,/) 1 / 2 . Our calculations are not precise enough to 
give critical exponents. 

In Fig. |l5], we show an approximate phase diagram in 
the U-h plane constructed using the calculations done 
on a single-size lattice (L = 30) for a single realization of 
the random potential with A/t = 1. We chose the first 
inflection point of the c^ -1 curve plotted against h as 
the rough phase boundary. As the interaction strength 
increases, the system becomes effectively less disordered, 
and therefore the critical value of h required to drive the 
system to the flux liquid phase becomes smaller. 

In Fig. [ll| we plot the curves for 9 V and c^ -1 versus 
h obtained for five different realizations of the random 
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potential. Due to the small system size (L = 30) used in 
the calculations, these quantities show sizable fluctuation 
effects. The relative fluctuation of the vortex tilt angle is 
seen to be smaller than that of the inverse tilt modulus. 

The Bose-glass to flux-liquid transition can also be 
driven by changing parameters such as the disorder 
strength A/t and the density n. In Fig. [l7|, we show 9 V 
and c v i4 ~ l as functions of A/t. A transition is observed 
at a critical strength of A/t. When the non-Hermitian 
field h is bigger, the transition occurs at a higher value of 
A/t. In this figure, we also compare the numerical results 
with the analytical formulae in the weak disorder limit 
derived in Appendix A. Since the variance of the ran- 
dom potential, (V 2 ) — (Vi) 2 , for the particular disorder 
configuration used in this calculation is approximately 
0.419, we relate the disorder parameter in the analytical 
formulae, T, to the parameter AbyTw 0.419A 2 . 

We have also computed the tilt angle and the tilt mod- 
ulus for the cases with a single attractive impurity. When 
both the interaction and the transverse external field is 
very weak, all flux lines can be pinned by one impu- 
rity and the system is an insulator. As the interaction 
strength is increased however, the screening effect causes 
the lowest-energy bound state to change into an extended 
one and the system becomes a tilted flux liquid. This is 
illustrated in Fig. ^ and the corresponding phase dia- 
gram is shown in Fig. 191 



V. CONCLUSION 

In the present paper, we have developed a numerical 
method of studying the strongly-disordered boson Hub- 
brad model in the presence of a constant imaginary vector 
potential and a random scalar potential and used it to 
understand the physics of vortices in high-temperature 
superconductors with columnar defects. We have found 
that the interaction causes the screening of a random 
potential, which is most effective for the ground and low- 
lying excited states. For sufficiently large values of the ef- 
fective interaction parameter or the non-Hermitian field, 
the nature of these states is changed from localized to 
extended. All physical properties are strongly influenced 
by these phenomena. In this paper, we presented a de- 
tailed study of the vortex tilt angle and the tilt modulus 
in (1 + 1) dimensions. 
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APPENDIX A: NON-HERMITIAN 
BOGOLIUBOV THEORY OF A DILUTE BOSE 
GAS IN A WEAK RANDOM POTENTIAL 

In the weak disorder limit, it is possible to obtain ana- 
lytical results for various physical quantities using the 
Bogoliubov approximation. In this Appendix, we de- 
scribe the results for both continuum and lattice models 
of a dilute Bose gas in a constant imaginary vector poten- 
tial and a weak random potential. In the section on the 
continuum model, we define V as the volume and n as the 
number of particles per volume. In the section on the lat- 
tice model, V represents the number of lattice sites and 
n is the number of particles per site. Our conclusions 
serve as a check on results obtained via coherent state 
path integrals in Ref. 7. The results for non-Hermitian 
lattice models with disorder are new. 



1. Continuum model 

We consider a dilute Bose gas of volume V in a weak 
random potential, V(x), and a constant non-Hermitian 
vector potential, g, in d dimensions. The Hamiltonian is 
given by 



d d x -ip^t/j^ipip 



' 2 



(Al) 



where ip(x) is the boson field operator and vq is the in- 
teraction strength. This model is a continuum version of 
the lattice model (|l|). Upon substituting 



V 



ik-x 

e a k , 



e lkx I4 



(A2) 



into Eq. (Al), we obtain 



x - (hk + igY t I ^ 

k k.q 



a k fl q 



v o sr^ t t 

2^ a k a k' a k"ak+k'-k" 



(A3) 



k,k',k" 



We assume that F(x) is a Gaussian random potential, 
characterized by a single parameter 7: 



y(x)y(y)= 7 J(x-y) i 

VkVq = 7<S k ,-q, 



(A4) 



where Vk=o = and the bar indicates averaging over 
disorder. 







Following Bogoliubov, we replace do and a\ by 

1 



'No 



N- 



(A5) 



and expand the Hamiltonian in terms of a k and a k with 
k ^ 0: 



n = v 



E 



2^+2 

h 2 k 2 



V n 



2m 



v n - 



%— k g) aja k 



k^O 

^E^K 

k^O 



(A6) 



where n = N/V. We ignore the cubic and quartic terms 
represented by ■ • • and diagonalize the remaining terms 
by the nonunitary canonical transformation 



2 1.2 



h z k 

2m 



h 2 k 2 

2m 



i—k ■ g 



i—k ■ g 

m 



(A7) 



where 



v n 



2l2 



2m 



'h 2 k 2 fh 2 k 2 



v n 



'n 2 k 2 fh 2 k 2 



2m V 2m 



2vqu 



2m \ 2m 

+ i — k ■ g. 

m 



2vqu 



The left and right annihilation operators c£ and sat- 
isfy bosonic commutation relations 

[<£, Cq] = 5 k , q , [c£, c«] = 0, [4, c$] = 0. (A9) 

The diagonalizcd Bogoliubov Hamiltonian has the simple 
form 



7~Lb = E q 



E, 



E WkC k C k' 
k^O 



9 9 , 1 2 

V 2to 2 



2 V ^ 



k/O 



1 E 



2 1.2 



w n 



'^ 2 fc 2 fh 2 k 2 



2m 

m 2 n 2 k 2 



2m V 2m 



V ^— ' aA-CL>_ k 2m ' 

k^O 



2uo?t. 
(A10) 



where B ff is the ground state energy and u; k is the Bo- 
goliu bov quasiparticle spectrum. From Eqs. (||), (Q) and 
( A1C| ), we derive the imaginary current, 



d d k 
(27r) d 



l^kl 



4 cos 2 



[n 2 k 2 /2m + 2v n + 2{g 2 /m) COS z I 
and the superfluid density, 



(All) 



n — 4n 



/d a k 



[h 2 k 2 /2m + 2v n + 2(g 2 /to) cos 2 0] ' 

8{g 2 /m) cos 2 
[Ti 2 k 2 /2m + 2v n + 2(g 2 /m) cos 2 0] ' 



(A12) 



where is the angle between k and g and 4> is the angle 
between k and Sg. In case of the Gau ssian random po- 
tential, the integrals in Eqs. (All) and ( A 12] ) can be per- 
formed easily. We summarize the results for the disorder- 
averaged imaginary current J/, expressed in terms of the 
dimcnsionless parameter g = g / y/mv^n, below: 



TO 



1-7 



1 



mvon 



1 



2v 2 n 2 \l h 2 (l+g 2 f/ 2 _ 

if d = l, 



>i = 



1-7 



2ttTl von g 



1 - 



1 



if d = 2 



(A13) 



(A14) 



The disorder- averaged normal fluid density n n (= n — n s ) 
is listed below: 



1 



2v 2 n 



mvon 1 — 2g 



7,2 



TO 



(l+3 2 ) 5/2 
1 + 2g 



if d = l. 



7,2 



1 



2nh 2 v o9 2 L(l + 5 2 ) 3/2 

if d = 2, <% || g, 

2 ( 1 \ 



(A15) 



(A16) 



n n = 7 



27r?i vq g 



7,2 



1 - 



if d = 2, <5g _L g. (A17) 



ferent method 



Eqs. (A15- Al7).have been derived previously using a dif- 



2. Lattice model 



We consider the d-dimensional non-Hcrmitian Hub- 
bard model for lattice bosons, Eq. ([!]), in a weak random 
potential V x . Upon substituting 
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ik-x 

e a k , 



k 



(A18) 



where V is the total number of sites and Vk = o = 0) m to 
Eq. (|), we get 



6r + iei — t cosh(g„a/?i) 



1 x t 
-j= Vk_ q a^a q 



a^ak 



i/ V tt 

^7 a k a k'ak"ak+k'-k", 



2V 



(A19) 



k.k'.k" 



where 



en = t^^cosh.(g u a/Ti) [1 — cos(k v a) 

V 

ej = t sinh(<7„a/?i) sin (fc„a) . 



(A20) 



Using the Bogoliubov prescription, we obtain an approx- 
imate Hamiltonian 



Ha = V -tn^2 cosh(g v a/K) + ^Un 2 
(en + iei + U n) a\.a^ 

k^O 

\ Un T. 



ak«-k + a ik a -k 



k#0 

-Vn ^2 Vk ( a k + a - k ) ' 

k^O 



(A21) 



where n is the number of particles per site. We diagonal- 
ize this Hamiltonian by the canonical transformation 



ak 



c R - a k c^ k 



u 



t _ c t - gkcgk 



where 



1 

Z7n 



[e R - iei) , 

k 

k (e R - iei) , 



(A22) 



+ Un - y/e R (e R + 2Un) 



The diagonalized Hamiltonian has the form 
H B =E g + Y^ ^k^c R , 



(A23) 



k^O 



7? 1 

-r^ = — tn^2 cosh(g v a/fi) + -Ur 



£ [ eR + Un ~ V^R (ch + 2J7n) 

k^O 

1 v-^ |Vk| 2 



V * — ' WkW-k 
k#0 



(A24) 



From now on, we will restrict our attention to the one- 
dimensional case. For a Gaussian random potential sat- 
isfying 



v x v y = r<5 x , y 



and in the thermodynamic limit, we obtain 



?±= S uxh(h)ll + ~ C — sin 2 (/. •/■_> 

Jo n . n 7T 



(A25) 



i r & 

2n 



CM + 2icosh(/i) sin 2 (fc/2) sin(fc/2) 



-r 



x dk 



tcosh(h) Un + tcosh(h) sin 2 (kj '2) 
cosh(2/i) + cos(fc) + 2 



2£/n cosh(/i) + i cosh(2/i) — t cos(fc) 



(A26) 



and 

= ncosh(h) — n n i — n„2 
"* dk ' 



n„i 



sin 2 (fc/2) cosh(/i) 



C/ncosh(/i) + 2icosh(2/i) sin 2 (fc/2) sin(fc/2) 



[i cosh(/i)] 1/2 f/n + t cosh(ft) sin 2 (fc/2) 



1/2 



1 



tsinh 2 (/i) [7n + 2tcosh(/i) sin 2 (fc/2) sin(fc/2) 



«n2 



= r 



[icosh(/i)] 3/2 C/n + i cosh(/i) sin 2 (fc/2) 
j£ f cosh(ft-) 3 cosh(2/i) + cos(fc) 



2£7 n cosh(/i) + i cosh(2/i) — t cos(fc) 



4 sinh 2 (/i) [Un + It cosh(/i)] cosh(2/i) + cos(fc) + 2 



2f7 n cosh(/i) + f cosh(2/i) — f cos(fc) 



(A27) 



where fc = fca. The quantity n n \ is the contribution 
to the disorder-averaged normal fluid density due to the 
breaking of the continuous translational symmetry in a 
lattice and n n i is the contribution due to disorder. 
After preforming the integrals, we find 

J i= sinHh) \ 1+ l 



ii 



1 

mr 

r 



■ tan 



cosh(/i) 



nir w cosh(ft) 

2 cosh(ft) + 2m + u sech 2 (ft) 1 
4 * 2 [u + cosh(ft)] 3/2 [u + sinh(/i) tanh(ft)] 3/2 J ' 

(A28) 



and 



n n i = —- cosh(ft) 



1 1 

H — 



\/ m cosh(ft) cosh(2ft) 



7T cosh 2 (ft) 
+ [cosh(ft) cosh(2/i) — usinh 2 (ft)] tan -1 



cosh(ft) 



1 tanh 2 (ft) J ■ 
7T u + cosh(ft) I 

+ [cosh 2 (ft,) - u 2 ] tan -1 



-u + cosh(ft) 



cosh(ft) 



»n2 

n 



t 2 \ 2 [u + cosh(ft)] 1/2 [u + sinh(ft) tanh(ft)] 1/2 
1 [2u cosh(ft) + cosh(2/i)] [u cosh(ft) + 2 cosh(2ft)] 
4 [u + cosh(ft)] 3/2 [u + sinh(ft) tanh(ft)] 3/2 cosh 2 (ft) 
1 [u + 2 cosh(ft)] [2m cosh(ft) + cosh(2/i)] tanh 2 (ft) 

^~ O 3/2 3/2 

^ [u + cosh(ft)] ' [u + sinh(/i) tanh(ft)] ' cosh(/i) 
1 [w + 2cosh(h)] 2 tanh 2 (/i) 

8 [u + cosh(/i)] 5/2 [u + sinh(/i) tanh(/i)] 5/2 cosh 2 (/i) 
x [8u 2 cosh 2 (ft) + 8ucosh(ft) cosh(2ft) 



- cosh(4ft) + 2] 



(A29) 



where u = Un/t. In the continuum limit where a — > 
0, ta 2 — > ft 2 /to a nd T a — > 7, th ese expressions reduce 
precisely to Eqs. (|A13| ) and (|a15| ). 
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Flux Liquid 




Bose Glass 



FIG. 3. Schematic phase diagram of high-temperature su- 
perconductors with columnar defects. H j_ is the transverse 
component of the external magnetic field. 



FIG. 1. One flux line (wavy curve) induced by the field H z 
and interacting with columnar defects in a cylindrical super- 
conducting shell. The trasnsverse component Hi is generated 
by the current I threading the ring. In this paper, we consider 
the situation where there are many interacting flux lines. 
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Re e x 

FIG. 2. Typical energy spectra of noninteracting bosons in 
a constant imaginary vector potential g and a random scalar 
potential in one dimension. 



FIG. 4. Occupation number N\ = (&a & a) of the Hartree 
eigenstates measured in the Bogoliubov ground state plot- 
ted against the winding number N w of the corresponding 
Hartree eigenfunctions. The parameters used are Un/t = 0.5, 
A/t = 1, h = 0.8 and L = 60. 
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0.08 




FIG. 5. Local density in the ground state no(i) for 
Un/t = 0.5, A/t = 1, h = 0.1 and L = 60. Inset: Comparison 
of the local ground state density with the random potential 
Vi. Vi is reduced by 30 times to make the comparion easier. 
Note that maxima in no(i) are anii-correlated with minima in 
the bare potential so that the screened potential is smoother 
than the bare one. 
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FIG. 6. Comparison of the screened random potential Wi 
with the bare random potential Vi for Un/t = 0.5, A/t = 1, 
h = 0.1 and L = 60. 
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FIG. 7. Autocorrelation function of the screened random 
potential for Un/t = 0.5, A/t = 1, h = 0.1 and L = 60 
compared to that of the bare random potential. 
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FIG. 8. (a) Hartree and (b) Bogoliubov spectra deep inside 
the superfluid (that is, flux liquid) region for Un/t = 0.5, 
A/t = l,h = 0.8, 0.6, 0.3 and L = 200. The same realization 
of the random potential was used for all plots. The linear 
slope of the Bogoliubov spectrum in the complex plane near 
the bottom of the band reflects the non-Hermitian analogy of 
the phonon part of a Bogoliubov spectrum. 
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FIG. 9. (a) Hartree and (b) Bogoliubov spectra near 
the Bose glass- superfluid phase boundary for Un/t — 0.05, 
A/t =l,h = 0.3, 0.44, 0.5 and L = 60. Plots for different val- 
ues of h are offset for clarity. There appears a small bubble 
of complex eigenvalues near the ground state in both spectra 
for h = 0.44. 
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FIG. 10. Participation ratios for the Hartree eigenstates 
plotted against the real parts of their energy eigenvalues for 
different system sizes. The parameters used are Un/t — 0.5, 
A/t = l,h = 0.8, 0.3 and L = 100, 200. 
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FIG. 11. Comparison of the Hartree spectrum with the 
winding number N w of the corresponding Hartree eigenfunc- 
tion. (a) Un/t = 0.5, A/t = 1, h = 0.3, L = 200. (b) 
Un/t = 0.05, A/t = \,h = 0.44, L = 60. In (b), the plot of 
the spectrum is offset for clarity. When the eigenfunction is 
real, we arbitrarily set N w = 0. 



FIG. 12. (a) Hartree spectrum and (b) screened potential 
in the presence of a single attractive impurity with Vd/t = — 5 
for Un/t = 0.07, 0.11, 0.15, h = 0.1 and L = 50. The impurity 
is located at i = 25. Plots of the spectra for different values of 
Un/t are offset for clarity. The lowest-energy extended state 
is chosen to have zero energy to align the bubbles. 
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FIG. 14. (a) Vortex tilt angle 9 V and (b) normalized in- 
verse tilt modulus C0/C44 (co = nei) for Un/t — 0, 0.015, 0.05, 
A/t = 1, n = 0.6 and L = 30 as functions of h. 
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FIG. 13. (a) Hartree spectrum and (b) screened potential 
in the presence of a single repulsive impurity with Vd/t — 5 
for Un/t = 0, 1.5,3,4.5, h = 0.1, and L = 50. The impurity 
is located at i — 25. Plots of the spectra for different values 
of Un/t are offset for clarity. 



17 



0.06 



0.04 



0.02 





Flux Liquid 
(Superfluid) 


Bose glass 









0.2 



0.4 



0.6 



0.8 



FIG. 15. Calculated phase diagram in the U-h plane for a 
single realization of the random potential with A/t = 1 on a 
30-site lattice, n (= 0.6) and t are fixed to constant values. 
Note that interactions shrink the region of Bose glass relative 
to the noninteracting case U = 0. 




FIG. 16. Fluctuations of (a) the vortex tilt angle and (b) 
the inverse tilt modulus for five different realizations of the 
random potential. The parameters used are Un/t — 0.05, 
A/t = 1, n = 0.6 and L = 30. 
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FIG. 17. (a) Vortex tilt angle and (b) inverse tilt modu- 
lus for Un/t = 0.025, h = 0.1,0.2, n = 0.6 and L = 30 as 
functions of A/t. Numerical results are compared with the 
analytical formulae in the weak disorder limit derived in Ap- 
pendix A. The disorder parameter in the analytical results, 
F, is related to the parameter A by T « 0.419A 2 . 
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FIG. 18. (a) Vortex tilt angle and (b) inverse tilt mod- 
ulus in the presence of a single attractive impurity with 
V d /t = -1.5 for Un/t = 0,0.005,0.015,0.05, n = 0.6 and 
L — 50 as functions of h. 
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FIG. 19. Calculated phase diagram in the U-h plane in the 
presence of a single attractive impurity with Vd/t — —1.5 on 
a 50-site lattice, n (= 0.6) and t are fixed to constant values. 



TABLE I. Correspondence of the parameters of the vortex 
line system with those of the boson system. 
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